Phase Transitions in Isolated Vortex Chains 



Matthew J. W. Dodgson 
Theory of Condensed Matter Group, Cavendish Laboratory, Cambridge, CB3 OHE, UK. 

(Dated: January 11, 2002) 

In very anisotropic layered superconductors (e.g. I^S^CaC^Orr) a tilted magnetic field can 
penetrate as two co-existing lattices of vortices parallel and perpendicular to the layers. At low 
out-of-plane fields the perpendicular vortices form a set of isolated vortex chains, which have re- 
cently been observed in detail with scanning Hall-probe measurements. We present calculations 
that show a very delicate stability of this isolated-chain state. As the vortex density increases along 
the chain there is a first-order transition to a buckled chain, and then the chain will expel vortices 
in a continuous transition to a composite-chain state. At low densities there is an instability to- 
wards clustering, due to a long-range attraction between the vortices on the chain, and at very low 
densities it becomes energetically favorable to form a tilted chain, which may explain the sudden 
disappearance of vortices along the chains seen in recent experiments. 

PACS numbers: 74.60.Ec, 74.60.Ge 



I. INTRODUCTION 

The vortex system in layered superconductors ,@ which 
includes the high-T c cuprates, displays a rich set of 
physical phenomena, such as a thermodynamic melting, 
transition^ a pinning-induccd disordering transitioiij 
and various structural transitions at different anglesD'E 
The average density and orientation of the vortices are 
controlled by the magnetic field, because each vortex car- 
ries one quantum of flux, <3>o = hc/2e. In this paper we 
are concerned with the vortex chains that appear in a cer- 
tain regime of tilted magnetic fields. These chains, which 
consist of a high density of flux lines perpendicular to 
the layers, were first observed with the Bitter decoration 
technique by -Bolle et al.B A qualitative explanation fol- 
lowed shortlyO in terms of the proposed crossing-lattice 
state in tilted fields.Utl This state consists of a lattice of 
flux lines perpendicular to the layers crossed by a lat- 
tice of flux lines along the layers. The in-plane lattice 
is strongly distorted due to the anisotropy,EB and has a 
large spacing along the layers. Huse surmised that a pos- 
sible attractive interaction between the two species of flux 
line would lead to a higher density of out-of-plane flux 
lines along chains, with an inter-chain separation equal 
to the in-plane flux-line spacing. This picture seems to 
be consistent with the experimental observations.^ The 
energy of the crossing-lattice state was considered by 
Benkraouda and LedvijJ13 who found a transition from 
a single lattice of tilted flux lines to the crossing-lattice 
state as the tilt angle is increased for sufficiently large 
anisotropy. Their work, however, neglected interactions 
between the two crossing lattices. 

Interest in the crossing-lattice state hasJieen revital- 
ized in the last couple of years. KoshelevEj has shown 
that the regime for crossing lattices is larger than previ- 
ously expected. This is because a correct treatment must 
include the interactions between the perpendicular flux 
lines, and a distortion away from ideal crossing lattices 
has a lower energy. Experimental evidence comes in the 
form of the unusual dependence of the melting tempera- 



ture of-the vortex lattice as a function of magnetic field 
angle,EilE£l which is consistent with the crossing-lattice 
state rather than a single tilted flux-line latticeO In 
addition, the work of Koshclcv explained quantitatively 
the attraction of perpendicular flux lines that leads to 
the vortex-chain state in a certain field range. Apart 
from more detailed work pp the melting transition^ and 
on magnetization curves these concepts have—also in- 
spired recent scanning Hall-probe measurementstS on the 
vortex-chain state. These experiments are similar to 
the Bitter decoration technique in that they only probe 
the field distribution emanating from the surface of the 
superconductor. However, they have the advantage of 
speed and resolution, and are not a "one-off" measure- 
ment, so that the system can be finely tuned to observe 
different effects. These new experiments have shown un- 
expected properties within the chain state. In particu- 
lar, it is possible to tune to a field range where all of the 
out-of-plane flux lines are arranged in chains. It is this 
"isolated-chain state" that we will study in this paper. 

We will use the following geometry. The z-axis is per- 
pendicular to the layers, and the component of mag- 
netic induction B z determines the density of "pancake" 
vorticesLH in every superconducting layer. The in-plane 
field is along the x-axis, and B x gives the density of the 
flux lines in this direction. These in-plane flux lines have 
their centers in the spacing between—layers, and so are 
commonly called Josephson vortices. i3 A z-directed flux 
line is made from a stack of pancake vortices, and con- 
tains circulating currents in the layers up to a distance of 
the penetration depth, A a b = A, from the vortex center. 
A Josephson vortex in the rr-direction has an ellipsoidal 
current pattern, with the flux and currents decaying over 
the much larger distance A c = 7A in the y-direction. 
Here, 7 is the anisotropy ratio, which is large for weakly 
coupled layers. The phase singularity of the Josephson 
vortex is confined between two neighboring layers, with 
separation d, where the phase difference across the layers 
changes by 2ir over a distance of the Josephson length 
jd. 
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FIG. 1: Optimal displacements of the vortex chain of pan- 
cake stacks due to the crossing of Josephson vortices, using 
the parameters appropriate for BiSCCO quoted in the text. 
The figure takes the results from calculations for a c h = 10A, 
and aj z = 20d which for the isolated-chain state corresponds 
to B x = 53 G and B z — 0.8 G. The displacements are magni- 
fied by two for clarity, and the z-scale is not the same as the 
a;-scale. 



The most simple structure of the crossing-lattice 
state is when an ideal triangular lattice (with spac- 
ing ap = [(2/v / 3)<l , o/-B z ] 1 / 2 ) of pancake- vortex stacks 
crosses a stretched triangular lattice of Josephson vor- 
tices, with separation in the z-direction of aj z = 
[(2/ V3)$o / Bx] 1 / 2 1 \fl ■ Koshelev has shown how the in- 
teractions between pancake stacks— and Josephson vor- 
tices lead to an effective attraction^ At the low pancake 
densities we are interested in here (B z < <I>o/A 2 ) this 
means that there will be large distortions from the ideal 
triangular pancake- vortex lattice, with a higher density of 
pancake stacks located over the centers of the Josephson 
vortices. This leads to the high-density chains ofpgapcake 
stacks that are observed by Bitter decoration,Btll with 
an inter-chain spacing of aj y — % /7[(v / 3/2)<i>o/B E ] 1 ^ 2 . 
Eventually, for small B z -C $o/A 2 , all of the pan- 
cakes lie over the Josephson vortices. This isolated-chain 
state has a pancake-stack separation along the chains of 
flch = (®o/B z )/aj y , see Fig.]!]. 

The crossing lattices must compete energetically with 
the more conventional tilted lattice of vortices. This is 
quite similar to the vortex lattice in a continuous super- 
conductor, only there is a kinked structure along each 
vortex with a periodicity L = d(B x /B z ). The anisotropy 
of the currents mean that there is a stretched aspect ra- 
tio of this lattice. In addition, this distortion raay be 
different from that expected from simple rescalingEa due 
to the attraction of tilted vortices along the direction of 
tilt. This may lead to a "tilted chain state" ,E3 which 
is distinct from the crossing chains mentioned above. 
Koshelev showed that for intermediate pancake densi- 
ties B z > $o/A 2 there should be a transition at a small 
value of B x = B* « 0.01B Z between the tilted lattice 
{B x < B*) and the crossing lattices (B* < B x ). 

In this work we will concentrate on very small out-of- 
plane fields and B z <C B x < $o/ A 2 when all of the pan- 
cake stacks become attached to the centers of the Joseph- 
son vortices, giving the isolated-chain state observed in 




FIG. 2: Phase diagram of the isolated vortex chain, show- 
ing four phases: Isolated chain of tilted stacks; Isolated chain 
of straight stacks crossing Josephson vortices; Isolated buck- 
led chain; state with chains plus ejected vortices. Axes show 
the spacing of pancake vortices a c h and the Josephson vortex 
spacing aj z within the chain. 




FIG. 3: Same phase diagram as in Fig. 2 but with axes 
converted to the the in-plane and out-of plane components of 
magnetic field, assuming the Josephson vortices form a regu- 
lar stretched triangular lattice. 



recent scanning Hall probe measurements. E3 We find that 
the stability of the isolated-chain state is quite delicate. 
We first note in Section |l| that the existence of a sta- 
ble crossing configuration in the isolated chains is only 
possible for large enough anisotropy, layer spacing, and 
pancake separation. We also find that the energy of the 
isolated chain as a function of pancake-stack separation 
has a minimum, i.e., there is an optimum density of pan- 
cake stacks along a chain, and there will be an instability 
towards clustering when the total density is low. As the 
density of pancake stacks along the chain increases (at a 
fixed density of Josephson vortices) , the chain will buckle 
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(Section |lll|), and then eject vortices at a critical mini- 
mum separation (Section |IV|). More surprisingly, at very 
low pancake density the chain state may be replaced by 
a chain of tilted stacks (Section [v]). These results are 
summarized in Figs. || and ^, where we plot the calcu- 
lated phase boundaries for three first-order transitions: 
clustering, buckling and tilting as well as a continuous 
ejection transition. Finally, in Section VI we discuss the 
extent to which these transitions have been observed ex- 
perimentally, and consider the effect of fluctuations due 
to finite temperature or quenched pinning disorder. 

We conclude this introduction with a note OB—pa^. 
rameters. The experiments we have referred toJjiiirLJ 
were performed on the extremely anisotropic cuprate 
Bi 2 Sr 2 CaCu20 I (BiSCCO). Where we make explicit cal- 
culations we will take the following appropriate param- 
eters 7 = 500, A = 2000 A, and d = 15 A. This gives 
a Josephson length (the size of the non-linear core of a 
Josephson vortex) of yd = 7500 A ^> A. 



II. STRUCTURE OF THE ISOLATED VORTEX 
CHAIN 

Koshelev has estimated the structure of a pancake 
stack that crosses the center of a single Josephson 
vortex£3 In the crossing configuration there is a Lorentz 
force on the pancakes in a direction parallel to the 
Josephson vortex, which must balance the attraction 
of each pancake to its stack. This distorts the pan- 
cake stack, with the largest displacements by the two 
pancakes immediately above and below the Josephson 
core. The distortions give an energy gain of AE X f» 
-8.4£ d(A/7d) 2 /ln(3.5 7 d/A). Here e d = ($ /47rA) 2 d is 
the typical energy scale for pancake interactions. This 
causes the pancake stack to be attracted to a stack 
of Josephson vortices within the crossing-lattice state. 
Koshelev's result uses a quadratic approximation for the 
pancake-to-stack attraction. Recent workEJ using the 
correct potential has shown that, while Koshelev's es- 
timate for the crossing energy is close to the full result, 
a stable crossing configuration only exists for extreme 
enough anisotropy, 7 > 7m i n = 2.86A/<i. In this section 
we include the interaction between pancakes in different 
stacks separated by a c h along an isolated chain. For fi- 
nite a c h we find a further reduction of the stability limit 
of crossing found in Ref. We also find a minimum in 
the energy of the chain as a function of stack separation, 
which leads to a clustering instability when a c h A. 

To calculate the crossing configuration of the chain, we 
use the following assumptions: First, we neglect the effect 
of induced Josephson currents due to displacing pancakes 
from their stacks (reasonable for large values oiyd/X and 
small enough displacements u < yd). Second, we utilize 
the long range of the remaining electromagnetic pancake 



interactions,!!!! of order A in the z-direction. There are 
many (~ X/d ~ 10 2 ) pancakes that contribute to the cur- 
rent distribution in one layer, determining the potential 



felt by a given pancake. We can therefore reduce this 
many-body optimization problem to a one-dimensional 
problem, considering the displacement of a row of pan- 
cakes in one layer under the potential due to the ideal 
chain of pancake stacks (the corrections due to the dis- 
placements in other layers is small when the number of 
pancakes with large displacements is much less than X/d). 

Within this scheme, the energy profile (per stack) for 
displacing by u n a row of pancakes in the nth layer when 
there is a Josephson vortex between layers and 1 cross- 
ing the chain is, 
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The first term here comes from the Lorentz force on 
the pancakes due to the in-plane current density from 
the Josephson vortex, J%, and tends to pull the pan- 
cakes away from their stacks. The form of this current 
will be discussed more in Section [II, but the numeri- 
cal value is taken from Ref. 



e.g., the current in the 
n = 1 layer (immediately above the Josephson vortex) 
is J\ — 2.28£oc/$o7^- The second term is the attrac- 
tive magnetic interaction of the pancake row with the 
remainder of all the stacks in the chain. Clem showed 
for a single pancake,Ea that this has the same interac- 
tion energy as the sum of a pancake with a full stack, 
V^ k (R) = 2e dK (R/X) plus a pancake with its anti- 
image VP^~ pml = 2eodln(R/L), so that for the entire row 
we find 



V:™(u) = 2e d[K (u/X)+]n(u/2X)+y E 



(2) 
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[Euler's constant y E is needed to fix the zero of energy to 
that of fully aligned stacks, c.f. the small x expansion of 
the modified Bessel function K (x)]. 

Such energy profiles for the n = 1 pancake row are 
shown in Fig. | for different values of a c h with the 
choice of the ratio yd/X = 3.75 (a reasonable choice for 
BiSCCO). For dilute chains there is a stable minimum 
(as expected from the results of Ref. ^j]), which deter- 
mines the optimal displacement. However, this stabil- 
ity disappears once the separation is below the critical 
value of a^ m = 5.3A. There can be no stable isolated-, 
chain state for densities higher than this critical value.t 2 ] 
Fig. |I] also suggests that the crossing configuration is only 
metastable, with a lower energy at ui — 0.5a c h, but this is 
not reliable as the simple Lorentz force argument does not 
hold for pancake displacements of such a large fraction of 
the stack-separation. For such displacements the original 
Josephson vortex becomes completely fragmented, and a 
different approach is needed. In fact, the competing state 
here is a "tilted chain" , and in Section [v| we will calcu- 
late the energy of the tilted chain, to compare to the 
isolated Qhain state. (A third possibility is a soliton-likc 
structures that we will not consider here.) 
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FIG. 4: Energy profiles for displacing a pancake row within 
a chain of pancake stacks, V^{ui), (dashed lines) and the 
total energy change AE\(ui) with a Josephson vortex be- 
low the pancake row (full lines), for different values of the 
pancake-stack spacing a c h. Note how the metastable mini- 
mum disappears for spacings smaller than a™ n = 5.3A. Each 
curve is only plotted up to iti = a c h/2. The fact that, for all 
values of o c h the lowest energy is for large displacements, is 
not a physical result, as the Lorentz force argument will break 
down when m is a significant fraction of a c h. These results 
are for a fixed value of jd/ X = 3.75. 
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FIG. 5: (a) The energy per pancake within the isolated 
chain as a function of separation o c h when aj z = 20d. Note 
the presence of a minimum energy at a c h = 8.2A. (b) The 
same result, but expressed as the energy density f c h{B z ) as a 
function of the out-of-plane flux density B z — $0 /a,chajy The 
thick line shows where the function has positive curvature, 
required to ensure local stability. The entire function must 
be convex, and so for densities less than B^ 1 the chain will 
phase separate into clusters with a c h = a *h an d regions with 
no stacks. 



The displacements of the entire stacks are shown in 
Fig. for a c h = 10A, and a Josephson-vortex separation 
of aj z = 20d (this corresponds to B x = 53 G for our 
parameters). All of the displacements contribute to the 
total crossing energy gain. We find that the crossing 
energy per pancake stack only lowers from AE x (oo) w 
-0.21e o rf for an isolated stack to AE X (a™ 11 ) w -0.26e o d 
as the stacks reach the critical separation a™ 11 . Therefore 
within the region of stable crossing, there is a pinning 
energy to the Josephson vortex centers per unit length of 
pancake stack, 

£jv- P in = AE x /a Jz 0.2e {d/aj z ), (3) 

where the last result is with our parameters for BiSCCO. 

In Fig. |^a we plot the total energy per pancake 
of the isolated chain E ch = \ £\ #0 V c s ^ ck (ja ch ) + 
{d/aj z )AE x (a c h). Note that there is a minimum at 
a c h = 8.2A, reflecting the fact that at large separations 
the pancake stacks attract each other. This unusual fea- 
ture has been explained by Buzdin and Baladiec3 by con- 
sidering the distorted pancake stacks of Fig. |l| as the su- 
perposition of a straight stack plus a series of pancake- 
anti-pancake dipoles. The straight stacks give a repulsive 
term, but this is exponentially small for a c h 3> A, while 
the dipoles give a weak attractive term that only falls 
off like (note there is some similarity to the at- 

traction between tilted flux linesc2l) . Therefore there is a 
net attraction at large distances, which will destabilize a 



homogeneous chain at very low densities. 

To formally describe this instability we should consider 
the free energy density as a function of out-of-planc field, 
f(B z ) = / C h(-Bz) + £stack5 z /$ - Here £ stack is the line 
energy of a pancake stack, f c h(B z ) = E ch (a ch ) / a ch aj y d 
and B z = $o/a cb aj y . In Fig. |b we plot f c h(B z ) using 
the same data in Fig. |^a. We also note that the thermo- 
dynamically stable phase is determined by the minimum 
Gibbs free energy g(H z ) = f[B z {H z )} - B z H z /Att with 
H z = Airdf /dB z . By a geometric construction we see 
that the point B z = B* 1 on f c h(B z ) where a straight line 
from the origin connects with the same gradient must 
have g = 0, and there is a first-order phase transition 
between a Meissner phase with no pancake stacks and a 
finite density of stacks with a c h = a*^. All the points 
at lower density B z < B* 1 have g > and so are ther- 
modynamically unstable. In fact, similar arguments de- 
termine that f c h(B z ) must always be a convex function. 
The dotted line represents the region where the curvature 
of f c h{B z ) is positive, meaning that the solution here is 
always unstable. 

The straight line that determines B* 1 is given by 
/ = BSH^/in where 8H*i is the change in the lower 
critical field due to the attraction of pancake stacks 
to the Josephson vortices. In Fig. |5|b it is given by 
5H* X = —2.5 • 10~ 3 H* 1 which is hardly measurable. In 
real experiments however, the geometry of the samples 
often makes demagnetization effects important such that 
the average B z becomes fixed to the external field H ext , 
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rather than to H. This means that small values of B z 
are accessible, and that for B z < B* 1 we may expect 
a coexistence of B z = and B z = B* 1 phases (c.f. the 
intermediate state in type-I superconductorsE3). Alter- 
natively we could describe this mixed regime in terms of 
"clusters" of pancake stacks with separation a*^. The 
relative proportion of space taken up by these clusters is 
determined by the value of B z , but the size of individual 
clusters depends on the energy of the "domain wall" be- 
tween the cluster of stacks and the region of no stacks, 
compared to the magnetic energy cost of large clusters 
with the wrong flux density. Also, in the experimen- 
tal situation one can have different pancake densities on 
the different chains, so the inhomogenous state may have 
coexisting empty and filled chains. Our results for the 
critical field separating the clustered phase from the ho- 
mogeneous isolated straight chains are shown in Figs. ^ 
and E3. 



III. BUCKLING INSTABILITY WITHIN 
VORTEX CHAIN 

It is important to realize that, while the isolated-chain 
state gains energy due to the Lorentz force of the Joseph- 
son vortex currents, there is an energy penalty for the 
pancake stacks to be close to each other. One conse- 
quence of this is a maximum density above which pan- 
cake stacks will be ejected from the chain. This criti- 
cal density will be derived in the next section. Below 
this density, however, the chain can already react to the 
stack repulsion by buckling. Note that we will assume 
that the Josephson vortices remain straight, as in exper- 
iments they are held in place by the interactions with 
Josephson vortices in neighboring chains, and the chain 
separation is much smaller than the interaction range for 
Josephson vortices, aj y <C 7A. In contrast, the interac- 
tions between pancake stacks on different chains are in 
the dilute limit, cij y S> A, and so are easily displaced. 

To calculate this buckling, we need to know the energy 
gain from a crossing event when the pancake stack is 
displaced away from the center of a Josephson vortex. 
We therefore need the full current profile of a Josephson 
vortex, which is a solution of the non-linear equations 
that arise from the London-Lawrence-Doniach model£rE£l 
An accurate numerical solution for a Josephson vortex is 
described in Appendix B of Ref. For a vortex directed 
along x. centered at y = and between the n =p_p, 1 layers, 
the current in the y direction can be writtenE3 (ignoring 
screening, i.e. n < X/d and y < A c ), 



(4) 



with p n (y) = 4>' n {y) the reduced superfluid momentum, 
where the phase has the form,Ej 
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FIG. 6: The energy profile when we pull one pancake stack 
away a distance u y from the isolated chain state, for a c h = 
5.5A, and aj z = 25d. The geometry is shown in the inset. 
The dashed line shows the energy from the crossing events on 
the displaced stack. The dotted line is the interaction energy 
of displaced stack with the other stacks on the chain. The full 
line is the sum of these two contributions. Note the presence 
of two minima, one for u y = 0, and a lower minimum at 
u y = 4.9A. 
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When y — 0, this gives the results used in Section O, 
with J n (0) = (2c£ /<Z>old)C n /\n - ||, and d = 0.57, 
C 2 = 0.86, C7 3 = 0.99, and C n « 1 for n > 3. 

We can now recalculate the crossing energy for an ar- 
bitrary distance between a pancake stack and the cen- 
ter of a Josephson vortex. Within the quadratic ap- 
proximation used in Refs.[l3] and ^8] one can easily solve 
the crossing configuration for a single stack. Writing 
AE n (u n ) — —(&od/c)J n (y)u n + \cm\, where we take 
a = (eod/X 2 ) ln(\/r w ) with r w a short distance cut-off, 
we find the relaxation energy, 



AE x (y) 



1 



ln(A/r„ 



£ [ Pn {yhd)f. (6) 



[in- 



i) 2 + y 2 + 0.38] 2 



Koshelev takes r w — u% (~ 0.29A 2 /7<i for y = 0), and in 
Ref. |2l| we have shown this to be a good approximation 
to the result when the full form of the pancake-stack 
interaction is used. 

This result for the crossing energy as a function of 
stack displacement u y is shown as the dashed line in 
Fig. 6, where there is one Josephson vortex for every 
25 pancakes along a stack. As might be expected, the 
energy increases quadratically at small u v . More inter- 
esting is the fact that the crossing energy vanishes only 
as AE X w —\ 2 /u y2 for u y ^d (note that an ex- 
ponential suppression will occur at extremely long dis- 
tances u v > A c ). In contrast, also shown in Fig. 6 
(the dotted line) for a c h = 5.5A is the repulsive energy 
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FIG. 7: (a) Energy of the buckled chain (see inset) as a func- 
tion of the buckle parameter u v , for stack separation along 
the chain a c h/A = 5.2, 5.4, 5.6, 5.8, and 6.0. For the separa- 
tion (Zch/A = 6.0 there is only one minimum at u v = 0. For 
fflch/A = 5.8 there is a second minimum at u y = 3.3A. The re- 
maining cases have their lowest minimum at finite u v . (b) The 
size of the optimum buckling distortion u" pt as a function of 

Ctch. 

cost for displacing one pancake stack away from a chain, 
£re P = E j¥ o Km ack (V0«ch) 2 + u^/a) . This contribu- 
tion vanishes exponentially when u y > VachA, and so the 
total energy is always negative for large enough u v . This 
gives the possibility of two minima in the total energy 
(see the full line in Fig. 6), and a first-order transition 
where the displacement jumps to a finite value as a c h de- 
creases. While we have started by considering displacing 
a single stack from the chain, it is this instability that 
drives a buckling of the whole chain. 

To estimate the buckling transition we use a variational 
approach for the buckled configuration (See the inset of 
Fig. 0a). This configuration is only characterized by a 
single displacement parameter it y , and so we can look at 
the energy profiles in u v for different a c h/A. First, the 
energy per pancake is calculated before crossing relax- 
ations, 




with yi = for i even, and y; = (—l)( t+1 " 2 u v for i 
odd. Next we checked for the stability of crossing, as 
in Section O, and found that the crossing configuration 
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FIG. 8: Energy density of the buckled-chain state with 
aj z = 20d as a function of out-of- plane field B z . Also shown is 
the result for the straight chain from Fig. 5b. The two points 
joined by the straight line have the same Gibbs free energy, 
and those densities between these points are thermodynami- 
cally unstable towards a mixed phase of straight and buckled 
chains. 



becomes unstable if a c h < 2.65A. Finally, to calculate 
the energy gain from crossing the Josephson vortices, we 
use the result (||) giving the total energy of the buckled 
chain, 

E h - c ( U y)=E*- c ( U y) + -?-[AE x (0)+AE x («»)]. (8) 

In Fig. 0a we plot this energy for a c h/A from 5.2 to 6.0, 
with aj z = 25d, showing how the minimum energy of a 
buckled chain crosses the energy of a straight chain as a c h 
decreases. Note that the minimum for a finite u v is quite 
shallow, and so we might expect large fluctuations in the 
extent of buckling along the chain due to random vortex 
pinning, which is discussed in Section yf\ . Also shown in 
Fig. 0b is the optimal displacement uj^as a function of 
a c h, which has a smallest value of «JJ pt = 2.3A at the first 
appearance of (meta)stable buckling at a c h = 5.8A and 
then increases at higher densities. 

It has been shown in Section || that one should plot 
the energy density of the chain as a function of the out- 
of-plane flux density B z , in order to determine the stable 
thermodynamic phases. Such a plot is shown in Fig. ^ 
comparing the straight and buckled chains. It shows 
that the isolated straight chains at B z = B* 2 ~ w 1.3 G 
has the same Gibbs free energy as the buckled-chain 
state at B z = B* z 2+ m 2.5 G. Therefore at fields with 
B* 2 ~ < B z < B* 2+ the thermodynamic phase will be 
a coexistence of straight and buckled regions, with the 
relative proportion linearly dependent on B z . Note that 
B* 2 ~ is only slightly higher than B* 1 so that there is 
only a small regime where pure isolated straight chains 
are the stable phase. The same procedure has been fol- 
lowed for a range of aj z /d, and the resulting boundaries 
to the mixed buckled chains are shown in Figs, pi and pi. 
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IV. EJECTION OF PANCAKE STACKS FROM 
ISOLATED CHAIN 



We now consider a simple model for the energy of the 
composite state where a fraction of the pancake stacks 
are not located on the chains. This model will show that 
there is a continuous phase transition between this com- 
posite state and the isolated-chain state as a function of 
B z . The order parameter of this transition is the density 
of the dilute lattice of pancake stacks not trapped on 
chains. The simple model assumes that the chain spac- 
ing ajj,, and the spacing between pancake stacks is much 
greater than the penetration depth A. In this case we 
can separate the total interaction energy density of pan- 
cake stacks in the form f tot = / c h(^ch) + /dii(^dii), where 
f c h and z^dii are the (2D) densities of pancake stacks on 
and off the chains respectively, and the total density is 
v c h + ^dii = ^tot = B z /<&q- The separation of pancakes 
along the chain is a c h = l/ajyi/^. The chain energy 
/ch(^ch) is to be calculated as in Sections || and III for 
straight and buckled chains respectively. We will see that 
the phase transition to the composite state occurs when 
there is a minimum in / c h(^ch), and so near this transi- 
tion we will expand, 



ichKh) = e'A 2 (^ch - ^ c c h) 2 



(9) 



The density of pancake stacks not on a chain is small 
near the transition, so that the interaction energy density 
is simply (using the limit Kq(x) w ^Jir/2xe~ x and only 
including nearest neighbors), 



/dil(^du) 



(10) 



Note that in the small i^dii limit, all terms in a power 
series expansion of /dn(^dii) are zero, i.e. the function is 
extremely flat. For this reason the critical total density 
at the transition is only determined by the minimum in 
/ch(^ch) at v^ h . For v tot < j/; h , all of the pancake stacks 
are on chains. At densities just above v^b the energy is 
minimized with a dilute off-chain density, 



fdu = 8v to 



3V27T eq 



1 



(11) 



The second term is extremely small for the fields we are 
interested in B z <C $ /A 2 , so that we can say in the 
composite state at z/ to t > there is a fixed density 
on the chains is c h = and the remaining density is 
The phase boundary r/L between the 



ftot 



composite state and the isolated-chain state is calculated 
from the minimum of / c h- In Figs. and ^| this bound- 
ary is plotted, but is not distinguishable by eye from 
the transition between mixed buckled and pure buckled 
chain. Therefore the pure buckled chain has a very nar- 
row range of existence. It is worth stating that au. ejec- 
tion transition has been observed experimentally^ with 
results that seem to be consistent with the above, al- 
though a quantitative analysis of the transition has not 
yet been published. 




FIG. 9: Configuration of the tilted chain, where pancakes in 
a stack are displaced from layer to layer by L — a c h/ajz- 



V. CROSSING-TO-TILTED TRANSITION 
WITHIN ISOLATED CHAINS 

In this section we calculate the energy of a chain of 
tilted vortices, and compare to the chain of crossing vor- 
tices. The tilted chain is an alternative configuration 
to the crossed chain with the same density of pancake 
and Josephson vortices (see Fig. ^) where the tilt angle 
is determined by tan# = L/d = <z c h/<2jz- These tilted 
stacks have an increased energy of pancake interactions 
but there may be an energy gain from replacing fully 
developed Josephson vortices with many short Joseph- 
son segments of length L < jd. To calculate the energy 
of the tilted chain, we use the linear papproximation of 
the London-Lawrence-Doniach modelpca where we ig- 
nore the non-linear effects from regions with large phase 
differences across neighboring layers. This approxima- 
tion is justified as long as the length L is much smaller 
than the Josephson length jd and the pancake separa- 
tion a c h- Note that the condition L < jd corresponds to 
B x /B z < 7, which is the case we are interested in. The 
energy due to pancake interactions in the linear approx- 
imation is, 



nrtilt 



1 f d 2 k dq 



[l + ( 7 A) 2 (fc 2 + g 2 )] \S Z 



8tt J (2tt) 2 2tt [1 + A 2 (fc 2 + q 2 )] [1 + ( 7 A) 2 fc 2 + AV] 

(12) 

with S z (k,q) = $odE 3 > e " I,T " le " !tl(j ° ch+, ' i) and <7 2 = 
(2/d) 2 sin 2 (gd/2). The integral should be cut off due to 
the pancake cores at k > 7r/£ a b and due to the layered 
structure at q > ir/d. Similarly, the Josephson segments 
interact with an energy contribution, 



TTitilt 



1 f d 2 k dq 



\S. r . 



8tt J (2ir) 2 2tt [1 + ( 7 A) 2 fc 2 + X 2 q 2 



(13) 



with 5^ ^sin(^)E 



; -ig(n+i)d e -ifc x [j'a ch + (n+i)L] 

The integrals can all be done exactly, leaving a sum 
over the values of k x = 27rm/a c h- We have evaluated 
these sums numerically to find the energy of the tilted 
chain. This is shown for aj z — 25d as the dotted line 
in Fig. [HI Also shown as the full line is the energy per 
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FIG. 10: Energy per pancake of the tilted chain (dotted line) 
and the chain of straight pancake stacks crossing Josephson 
vortices, for aj z = 25d, plotted against the pancake separation 
a c h/A. Also shown (dashed line) is the difference of the two 
energies, which crosses zero at a c h = 45 A. 



pancake of the crossing chain with the same density of 
pancakes and Josephson vortices, Serosa = ^ch + Ef v oss , 
where E c ^ is defined in Section O and 



across 



Sod a c h 
7 ajz 



1.55 



(14) 

The tilted chain has a lower energy from Josephson 
vortex contributions (dominant at low density) but a 
higher energy from pancake vortex contributions (dom- 
inant at high density), and there is a first-order tran- 
sition between the isolated chain of straight stacks and 
the chain of tilted stacks. In fig. [l(] it is shown that 
for a c h < 45A (higher density) the preferred state is the 
crossing chain, while for a c h > 45 A (lower density), the 
lowest-energy is for the tilted chain. Note that Savelev 
et al.J^J also find a "reentrant" transition back to a tilted 
lattice for fields close to the a6-plane, when A > 'yd. 

A strict treatment should follow Sections || and III and 
consider the points of the same Gibbs free energy den- 
sity. This is not possible within our level of treatment, 
however, as we would need accurate results at very low 
B z < B x /j where non-linear effects in the tilted chain 
are important. Here, we can only show that there must 
be a low-density transition to the tilted chain, but we 
cannot plot the coexistence boundaries. Even so, we plot 
the line where the energies cross in Figs. || and || sepa- 
rating the clustered chain from a tilted chain phase. In 
real experiments there is a finite 2D density of Josephson 
vortices, which cannot be considered isolated in the sense 
of the pancake vortex chains (the interaction range A c is 
greater than the chain spacing aj y ). Therefore a full cal- 
culation of the 3D tilted lattice requires a minimization 
over the aspect ratio, which we leave for a future work. 



VI. DISCUSSION 

We now discuss the extent to which these transitions 
have been observed experimentally. In the receut. scan- 
ning Hall probe experiments of Grigorenko et al.,E£l there 
does seem to be a sharp transition as a function of B z 
between a composite state of chains with a dilute lattice 
and the isolated chain state, as derived in Section IV. 



At lower fields there is some evidence for buckling and 
clustering (Sections [II] and [II) although the influence 
of pinning disorder may be contributing to this. Finally, 
Grigorenko et al. report a strange transition at very low 
B z where the chains are replaced by faint, homogeneous 
"stripes" of flux. These stripes would be consistent with 
the flux distribution from a tilted chain, and it is likely 
that a transition from the isolated-chain state to a tilted- 
chain state (as found in Section [y]) takes place. Experi- 
mentally, a large jump in B z (i.e. in the pancake density) 
is seen at the transition to tilted chains; unfortunately we 
were not able to calculate the jump in this paper, as it 
requires an accurate treatment of non-linear effects in the 
small B z limit. Other experimental features are that the 
tilted chains are much straighter than the crossing chains, 
and less pinned. This is consistent with the fact that 
the tilted chains are really part of a 3D line lattice with 
stronger interactions that dominate over pinning disor- 
der, whereas in the crossing case the pancake stacks are 
in the extreme dilute limit, and therefore easily pinned 
in random low-energy sites. 

Finally, we briefly discuss the effect of fluctuations. 
The calculations in this paper have found the lowest- 
energy vortex configurations of the chains. In reality 
there will be distortions of these states due to thermal 
fluctuations and random pinning to inhomogeneities in 
the underlying crystal. It is well established that ther- 
mal fluctuations can lead to a melting of the vortex 
crystalll A melting transition to completely decoupled 
layers is studied for B x = in Rcf. BlJ involving short- 
wavelength fluctuations of pancake stacks with k z > n/X. 
While there will be some modification in the presence of 
Josephson vortices in tilted fields, in the low out-of-plane 
fields studied here, B z <C &o/\ 2 , this melting will take 
place at a high temperature —close to the Berezinskii- 
Kosterlitz-Thouless transitions at Tbkt = eod/2. If we 
consider long wavelength fluctuations, k z — > 0, then we 
may expect some kind of melting of the chain at very low 
fields due to the exponentially small interaction of pan- 
cake stacks. However, the mechanism of such a transition 
must be different from the "entanglement" proposed in 
the original-derivation of low-field vortex lattice melting 
by Nelsonrj Considering the nature of the crystalline 
order of the isolated vortex chain, we should recognize 
the two-dimensional nature of this state. A simple con- 
sequence is that thermal fluctuations at long wavelength 
will lead to a quasi-long-range ordered state. It may then 
be possible to have a 2D continuous melting transition via 
the unbinding of dislocations. However, a dislocation for 
this system of stacks corresponds to a stack that termi- 
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nates in the middle of the chain, and this should cost an 
energy linear in the system size, rather than the usual 
logarithmic energy of a dislocation. This will suppress a 
dislocation unbinding transition. 

In the experiments of Grigorenko et al. the pinning 
disorder seems to be more important in disturbing the 
chain states than thermal fluctuations. The fact that the 
pancake stacks are observed in fixed positions is due to 
pinning (otherwise thermal fluctuations would smear out 
the average density in the chain state), and this also tends 
to disorder the chains. On general grounds we expect 
pinning-induced wandering of the pancake stacks within 



the chains so that there is only a short range order. There 
will also be significant transverse displacements, and the 
buckling effect should be enhanced by disorder. 
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Ziirich was financially supported by the Swiss National 
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